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Abstract. We analyze longitudinal pion spectra from y^ijvjv = 2 GeV to ^/snn = 20 GeV within Landau's 
hydrodynamical model and the UrQMD hybrid approach. From the measured data on the widths of the 
pion rapidity spectra, we extract the sound velocity c 2 . in the dense stage of the reactions. It is found 
that the sound velocity has a local minimum (indicating a softest point in the equation of state, EoS) 
at y^ijvjv = 4-9 GeV, an energy range accessible at the Facility for Antiproton and Ion Research (FAIR) 
as well as the RHIC-Beam Energy Scan (RHIC-BES). This softening of the EoS is compatible with the 
formation of a QGP at the onset of deconfinement. 

The extracted sound velocities are then used to calculate an excitation function for the mean transverse 
mass of pions from the hybrid model. We find that, above -^/sjvjv ~ 10 GeV, even the lowest c? a gives a 
considerably larger < mr > of pions compared to data. 
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Over the last years, a wealth of detailed data in the 
Eiab = 20A — 160A GeV energy regime has become avail- 
able. The systematic study of these data revealed surpris- 
ing (non-monotonous) structures in various observables 
around Ei^ = 30A GeV beam energy. Most notable ir- 
regular structures in that energy regime include pQ, 

— the sharp maximum in the K + /ir + ratio [HEJH], 

— a step in the transverse momentum excitation function 
(as seen through (mi) — mo ) [3JS1II], 

— an apparent change in the pion per participant ratio 

® 

— a sign change in the v\ slope [B] 

— an apparent maximum in the final state eccentricity 
from an Hanbury Brown Twiss (HBT) analysis [7] and 

— a minimum in the moment products of net-proton dis- 
tributions (<t 2 k) [8] 

It has been speculated, that these observation hint 
towards the onset of deconfinement already at energies 
around y/s^N = 5-10 GeV. Indeed, increased strangeness 
production [S] and a non-monotonous fluctuation signal 
have long been predicted as a sign of a mixed phase at 
the onset of Quark Gluon Plasma (QGP) formation [TU1 
rTTUT5HT3 II 1 41115] within different frameworks and observ- 
ables. The suggestion of an enhanced strangeness to en- 
tropy ratio (~ K/tt) as indicator for the onset of QGP 
formation was especially advocated in [IB]. Also the high 
and approximately constant K inverse slopes of the tut 



spectra above ^/snn ~ 7 GeV - the 'step' - was also 
found to be consistent with the assumption of a parton 
hadron phase transition at low SPS energies [T71ll8j . 
Surprisingly, transport simulations (supplemented by lat- 
tice QCD (1QCD) calculations) have also suggested that 
partonic degrees of freedom might already lead to visible 
effects at ~ 7 GeV [HUMI] ■ Finally, the comparison of 
the thermodynamic parameters T and [Ib extracted from 
the transport models in the central overlap region [22] 
with the experimental systematics on chemical freeze-out 
configurations [23,24,25] in the T — fiB plane do also sug- 
gest that a first glimpse on a deconfined state might be 
possible around ^snn = 5-10 GeV. 



0.1 Longitudinal Dynamics 

In this paper, we explore whether similar irregularities are 
also present in the excitation function of longitudinal ob- 
servables, namely rapidity distributions. Here we will em- 
ploy Landau's hydrodynamical model p51l271B51BgpH l3lH 
[32] and the more sophisticated UrQMD hybrid approach. 
Landau's model entered the focus again after the most 
remarkable observation that the rapidity distributions of 
newly produced hadrons at all investigated energies can 
be well described by a single Gaussian at each energy. 
The energy dependence of the rapidity width can also be 
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Fig. 1. The root mean square a y of the rapidity distribu- 
tions of negatively charged pions in central Pb+Pb (Au+Au) 
reactions as a function of the beam rapidity y p . The green 
short-dashed line indicates the Landau model prediction with 
Cg = 1/3, while the grey lines with symbols indicate the 
UrQMD hybrid model results with different equations of state. 
Data (open symbols) are taken from [5ll48] . The statistical er- 
rors given by the experiments are smaller than the symbol 
sizes. Systematic errors are not available. 



reasonably described by Landaus's model. For recent ap- 
plications of Landau's model to relativistic hadron-hadron 
and nucleus-nucleus interactions the reader is referred to 
[55U54TI551I551I37] (and Refs. therein). Due to the simplicity 
of the Landau model we treat it as an approach to extract 
the qualitative features of the result. 

In several publications [381139] it has been shown that 
the results of Landau's model are not unambiguous and 
more sophisticated models yield quite different relations 
between the longitudinal expansion and the effective equa- 
tion of state. Therefore, to complete the investigation we 
will also employ the UrQMD hybrid model which has been 
developed recently to combine the advantages of hadronic 
transport theory and ideal fluid dynamics |40) . It uses ini- 
tial conditions, generated with the UrQMD event gener- 
ator pLTlFL^] . for a full 3+1D ideal fluid dynamical evo- 
lution, including the explicit propagation of the baryon 
current. After a transition back to the transport descrip- 
tion, the freeze out of the system is treated dynamically 
within the UrQMD approach. The hybrid model has been 
successfully applied to describe particle yields and trans- 
verse dynamics from AGS to LHC energies |4"01l4"51l44ll451 
[46] and is therefore more reliable than the basic approach 
by Landau. One major advantage of the hybrid model is 
that, in contrary to Landaus's model, the transparency of 
the baryon currents is fully taken into account as done in 
the UrQMD pure transport calculations eliminating this 
uncertainty. 

The main physics assumptions of Landau's picture are 
as follows: The collision of two Lorentz-contracted nuclei 
leads to full thcrmalization in a volume of size V/'jcMi 
"/cm being the Lorentz factor in the center of mass frame. 
This justifies the use of thermodynamics and establishes 



the system size and energy dependence. For simplicity, 
chemical potentials are not taken into account in this ap- 
proach. From these assumptions follows a universal for- 
mula for the distribution of the produced entropy, de- 
termined mainly by the initial Lorentz contraction and 
a Gaussian rapidity spectrum for newly produced parti- 
cles. Further decays of resonances are not included in this 
approach. Under the condition that the speed of sound, 
c s is independent of temperature, the rapidity density is 
given by [29 j f30 ] : 
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where if is a normalization factor, m p is the proton mass 
and ^snn the center of mass beam energy. Under these 
assumptions, Landaus's model allows to relate the ob- 
served particle multiplicity and distribution in a simple 
and direct analytic way to the parameters of the QCD 
matter, namely the equation of state, under considera- 
tion. 

For the hybrid model such a connection cannot easily be 
made and the width of the rapidity distributions may have 
non-trivial origins. In particular the initial state plays a 
role, as one assumes an early system which may not be 
in equilibrium. Also the length and therefore importance 
of the fluid phase becomes larger with increasing energy. 
Consequently we will perform several calculations with 
the hybrid model, fixing the EoS to different values of 
the sound velocity, in the fluid dynamical phase. Then 
we extract the mean root square width numerically from 
the final state particle properties and extract the needed 
sound velocity from a comparison to experimental data. 

As the hybrid model is a rather complex approach, it 
contains several parameters which can have an influence 
on the resulting values of the rapidity widths, e.g. the 
hydro-start time or the transition criterion back to the 
hadronic cascade. Studies on the parameters have been 
performed in [40 and have have shown a rather small 
influence on results. However, further and more detailed 
studies, applying statistical methods which exceed the scope 
of this paper, are in order and are currently under way 
(47j . One effect which might have a considerable impact 
on our results, the contribution of resonance decays on 
the pion spectra, will be discussed in the work. We will 
estimate the resonance contribution in the hybrid model 
and compare with the effect of the equation of state. 

Let us now analyze the available experimental data on 
rapidity distributions of negatively charged pions in terms 
of the hydrodynamical models. Fig.Q]shows the measured 
root mean square a y of the rapidity distribution of nega- 
tively charged pions in central Pb+Pb (Au+Au) reactions 
as a function of the beam rapidity in the center of mass 
frame. The green short-dashed line indicates the Landau 
model predictions with the commonly used constant sound 
velocity c 2 s = 1/3. The grey lines with symbols depict the 
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Fig. 2. Speed of sound as a function of ^/sjvjv for central 
Pb+Pb (Au+Au) reactions as extracted from the data using 
Eq. (open symbols) and from comparing to hybrid model 
results (full symbols). The statistical errors (not shown) are 
smaller than 3%. 



UrQMD hybrid model results with different equations of 
state (different but constant values of c 2 ), while the data 
points [51H8] are depicted by open symbols. 

Note that the black solid line in figure [T] depicts re- 
sults where the hybrid model is applied using an equation 
of state of an hadron resonance gas. When we suppress 
the resonance contributions to the pion spectra we obtain 
the dashed black line. From this comparison we can con- 
clude that resonance decays only play a minor role in the 
width of the rapidity spectra. Even more, the effect is only 
a shift to larger widths and does not change our results 
qualitatively. 

At a first glance the energy dependence looks struc- 
tureless. The data seem to follow a linear dependence on 
the beam rapidity y p without any irregularities. However, 
the general trend of the rapidity widths is also well repro- 
duced by the hydrodynamic model. However, there seem 
to be systematic deviations from an energy independent 
speed of sound. At low AGS energies, the experimental 
points are generally under predicted by Eq. ([2]), while in 
the SPS energy regime Landau's model overpredicts the 
widths of the rapidity distributions. Exactly these devia- 
tions do allow to gain information on the equation of state 
of the matter produced in the early stage of the reaction. 
For Landaus's model one can analytically invert Eq. ^ 
and express the speed of sound c 2 in the medium as a 
function of the measured width of the rapidity distribu- 
tion: 





(3) 

From the hybrid model we can extract the average 
speed of sound by a comparison of the experimental re- 
sults with the width obtained from different equations of 
state. 




\, [GeV]° 

Fig. 3. Excitation functions of the mean transverse mass for 
different fixed values of the speed of sound (grey lines with sym- 
bols). The black symbols denote the mean transverse masses 
with speeds of sound that were extracted from the rapidity 
distributions in the previous section. The black lines corre- 
spond to hybrid model results with a hadron resonance gas 
EoS. Here we compare cases where final interactions are sup- 
pressed (dashed line) and contributions from resonance decays 
were omitted (short dashed line). Data are depicted as open 
triangle and square symbols. 



Let us now investigate the energy dependence of the 
sound velocities extracted from the data. Figure [5] shows 
the speed of sound as a function of beam energy for cen- 
tral Pb+Pb (Au+Au) reactions as obtained from the data 
using Eq. ([3]) and Landau's model (open symbols) as well 
the speed of sound which was extracted from the full 3+1 
dimensional hybrid model simulations (full symbols). One 
can clearly see that the absolute values of c 2 depend con- 
siderably on the framework that is used. In general the 
speed of sound extracted from the full 3+1D hydrody- 
namics is smaller than from Landau's model. On the other 
hand, the sound velocities for both approaches exhibit a 
minimum (usually called the softest point). For the Lan- 
dau case it is not very pronounced and around a beam 
energy of ^/snn = 6-9 GeV while for the hybrid model 
it is more pronounced and as at -Jsnn — 4 GeV. The 
observed shift in the softest point can be explained by the 
following considerations. In Landau's model the full evolu- 
tion is governed by the hydrodynamic flow, which starts to 
develop once the two colliding nuclei have overlapped. In 
particular there is no pre-equilibrium dynamics. In the hy- 
brid model we observe a dynamical equilibration process 
and strong pre-equilibrium dynamics. This is particularly 
important at low energies where the length and relevance 
of the fluid dynamical phase is the smallest, when com- 
pared to the initial and final decoupling phase. The ex- 
tracted value of c 2 , of the hybrid model, is therefore only 
sensitive on the most dense part of the collision. Therefore 
one should regard the two results as bounds on the energy 
range in which the softening of the EoS occurs. 
A localized softening of the equation of state is a long pre- 
dicted signal for the formation of a Quark Gluon plasma. 
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Possible scenarios for for the observed softening include 
the presence of a mixed phase or a supercooled QGP in 
the specified energy range 49 ,50,51,52, 53, 54j . Therefore, 
we conclude that the measured data on the rapidity widths 
of negatively charged pions are indeed compatible with the 
assumption of the onset of deconfinement in the FAIR en- 
ergy range. However, presently we can not rule out that 
also an increased resonance contribution may be the cause 
of the softening [55] . 

0.2 Transverse Dynamics 

For a meaningful interpretation of the extracted speed of 
sound we relate the results found in the previous section to 
bulk transverse observables, in particular the mean trans- 
verse mass of pions. Fig [3] displays the excitation functions 
of the mean transverse mass of charged pions from the 
hybrid model with different fixed values of the speed of 
sound (grey lines with symbols), compared to data (open 
symbols). The black open circles depict mean transverse 
masses one would obtain with speed of sound extracted in 
the previous section. 

From this figure it is clear that in the AGS energy 
range the mean transverse mass is almost independent of 
the applied speed of sound, while at higher energies even 
a very soft equations of state does not allow to obtain a 
sufficiently small < mr > as has been measured. To un- 
derstand this behavior, figure |3] also shows hybrid model 
results where we applied an EoS that resembles a hadron 
resonance gas (black lines). Here we compare results of 
the standard implementation of the hybrid model (black 
solid line) with results where the final cascade stage is 
omitted (black dashed line) and the feed down contribu- 
tion of resonances to the pion mean transverse mass is 
omitted (black dotted line). One observes that the final 
cascade stage does not have any considerable influence on 
the mean < mj- >, while the resonance feed down con- 
tribution is of the same order than the dependence on 
the equation of state. In this context the mean transverse 
mass of pions seems not a well suited observable to pin 
down the equation of state, as one does not have full cer- 
tainty on resonance contributions to the results. Therefore 
we advocate that the rapidity spectra might allow for a 
clearer investigation of the EoS due to the smaller effect 
of the resonance decays. 

0.3 Summary 

In conclusion, we have explored the excitation functions of 
the rapidity widths of negatively charged pions in Pb+Pb 
(Au+Au) collisions. 

— The rapidity spectra of pions produced in central nucleus- 
nucleus reactions at all investigated energies can be 
well described by single Gaussians. 

— The energy dependence of the width of the pion rapid- 
ity distribution follows the prediction of hydrodynami- 
cal models if a variation of the sound velocity with the 
energy is taken into account. 



— The speed of sound excitation function extracted from 
the data has a pronounced minimum (softest point). 
Although the energy at which we extracted the softest 
point depends on the hydrodynamical framework, its 
existence turns out to be rather robust and it is lo- 
cated in the energy range of i?i a b = 5A-30A GeV. 

— This softest point is compatible with the formation of 
a Quark Gluon Plasma indicating the onset of decon- 
finement at this energy. 

— Transverse observables like the mean transverse mass 
of pions are shown to strongly depend on the contri- 
butions of resonance decays and are therefore a more 
ambiguous tool to discriminate different equations of 
state. 

Further explorations of this energy domain is needed and 
can be done at the FAIR and by CERN-SPS and BNL- 
RHIC experiments. 
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